Application of the Lorentz-Transform Technique to Meson Photoproduction 
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We show that the Lorentz integral transform ( LIT ) technique which has been successfully applied 
' ' to photoreactions in light nuclei can also be applied to photoreactions involving particle production. 

f--^ A simple model where results are easily calculable in the traditional fashion is used to test the 

f~>) , technique. Specifically we compute inclusive tv + photoproduction from deuterium for photon ener- 

^s^j ■ gies less than 200 MeV using a Yamaguchi model for the NN interaction. It is demonstrated that 

although the response functions for inclusive meson production do not have favourable asymptotic 
behavior one can nontheless extract them by inversion of the transform. The implication is that 
J*r . one can treat realistic problems of photo-meson production including all final state interactions by 

' means of the LIT technique. 

t> : 

PACS numbers: pacs numbers: 13.601e,21.45.+v,25.10.+s,25.20.1j 

<N ■ 

> : 

W> , I. INTRODUCTION 

in 

^vq ■ In a series of papers 0, IE IE EL IE IE it has been demonstrated that the LIT technique allows a convenient 
calculation of inclusive photoreaction cross sections wherein final state interactions are fully included. Further, the 
technique has been extended [EEI to exclusive reactions. A merit of the technique is that the calculation of continuum 
wave functions is avoided. In fact the differential equations to be solved are inhomogeneous and have solutions bounded 
, at infinity. The work cited above is based on non-relativistic quantum mechanics and a nucleons-only subspace. 

In order to treat meson photoproduction from nuclei it is desirable to see if the LIT technique can be extended to a 
larger subspace i.e. one including nucleons and mesons. This would enable the full inclusion of both meson-nucleon and 
NN interactions in the production calculations. As a test we consider inclusive photoproduction of low energy ( < 40 
MeV ) 7r + mesons from deuterium. A traditional calculation of this process was reported by Dressier, MacDonald, and 
Q'Connell[Tol| and by Gupta, Anand, and Bhasin [nj. Both of these groups employed a simple Yamaguchi potential 
|12| with parameters chosen to account for deuteron binding and low energy NN scattering properties. We adopt this 
model with the parameters used in 01 • Thus this example attempts to apply the LIT technique in the NN7T subspace 
where for simplicity only an NN interaction is included. Although no 7r— N scattering is contained in this model, this 
is not essential for testing the method. Our results will show that it is indeed possible to apply the LIT method to 
this problem and to extract the total photoproduction cross section. We should point out that our restriction to a 
simple NNpotential for this test calculation does not modify these conclusions. This is clearly seen from the series of 
papers 0, d IE EL IE IE where a range of NN potentials, from simple phenomenological models to modern realistic 
potentials, have been used in 3-nucleon and 4-nucleon problems. In no case did the complexity of the NN potential 
model have any effect on the accuracy or implementation of the LIT technique. With assurance that the formalism 
works the addition of a meson-nucleon interaction can be viewed as a technical point, the consequence of which would 
be more complicated numerics. 

Low energy photoproduction is described by the Kroll-Ruderman operator 

n int {x) = -ie(^\y^e x -a<J)r-{j)e- ikx <) )+ {x)5{x-x j ) . (1) 
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where k and ex denote the incident photon momentum and polarization vectors , respectively, and 4>+(x) denotes the 
meson field operator 

' e iq - x aUq) + e- iq - x a-{q)] (2) 
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with qo the energy of the meson and and a are the usual creation and annihilation operators. Here e is the positive 
elementary charge, is the 7r + -meson mass, ,f is the 7T-N coupling constant, a denotes the Pauli spin matrix vector 
and r_ the isospin operator. Dressier et al. [jjj] show that the KR term gives nearly the entire cross section for pion 
energies in the range considered here. Due to cancellations other terms in the production operator i.e. terms required 
by gauge invariance, do not contribute significantly at low pion energies. With E cm as the incident cm. energy the 
response function lZ(E cm ) for this process is 

Tl(E cni ) = 1 ££|UI 6(k, A) | D, M d ) | 2 5{E cm - E f ) , (3) 

M d ,\ f 

where | D, Md ) denotes the deuteron ground state with polarization Md and 

6(k,\) = Jd 3 qO(k,X,q)al(q) (4) 

with 

0(k,X, ql = -ie (X) _J= ^^^,(,VVN . (5) 

In Eq. (PJ above | /) denotes the wave function of the relative motion of the nnir + system with energy Ef, and 
indicates an integration over all relative momenta and a sum over all final nucleon spins. The inclusive cross section 
is related to TZ(E cm ) by 

2tt 2 

<y(E C m) = —n(E cni ) . (6) 

Relative momenta in the final state are taken as 

Px = (pi-p 2 )/2 (7) 

m * f r \ i 2m n , Q S 

Py = ^^4^ Vl+P2 > + ^f p ^ 

where p\,p2, and are the momenta of the final state neutrons and pion, respectively, and M is the total mass 
2m n + m-x with the mass of the neutron m n . As implied by this separation we are restricting ourselves to non- 
relativistic kinematics for both the nucleons and the pion. This allows the separation of the hamiltonian into Jacobi 
coordinates for the NN-7T three body system. In terms of these quantities the energy conserving 5-function appearing 
in TZ(E cm ) takes the detailed form 

5(E cm -Ef)=S\E CID -(2m n + m n )-^-^-] , (9) 

where fi is the reduced two-neutron-pion mass. Finally it is more convenient to use the cm. energy above threshold 

k 2 

W = k+ (2m„ + m 7! - m d ) > 

2m d 

where md is the mass of the deuteron. The LIT, referred to hereafter as the transform, of the response function 1Z(W) 
is then defined with an, oj > as 



L{<7R,cri) = / 
Jo 



'dw R W , (10) 

(W - <7 R ) 2 + erf 



dW < D, M d | 6\k, A) J {W g) 2 6(k, A) | D, M d ). (11) 



Because of the relation of k to W, a straightforward integration of the above equation would leave the operator 0(k, A) 
depending in a complicated way on the Hamiltonian H. Therefore, we proceed by setting k appearing in the operator 
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O to a constant arbitrarily chosen " pseudo-momentum" fc p . As a result we introduce a new transform Lfc p , which 
depends on k p and takes the form 



L K {a R ,ai) = \ ]T (D,M d \6\k p ,\) 1 6(k p , A) | D, M d ] 

M d ,X ^ 17 R * 111 



M d ,\ 

where the Lorentz function ip is solution of the inhomogeneous equation 

(H-a R + io!) | i> Md ,x((^R, <ti, fcp) > = 0(k P , A) | D, M d ) 

with 



[p x ,Py\H\pJ,p') 



pI , p 2 



KPy-Py) 



(12) 
(13) 

(14) 
(15) 



The inverse to the transform Li~ p (gr, ui) will lead to a corresponding response function TZ^ (W) and in turn will 
yield the correct cross section for k p = k. Consequently for each photon energy the calculation of the transform and 
its inversion must be repeated. 

As mentioned earlier the simple model here only includes an NN potential, i.e., V(p x ,pJ). The kinetic energy of 
the meson with respect to the nn pair appears in H as the term proportional to p y . Meson rescattering could be 
included by adding to H appropriate potentials. Although that would considerably complicate the numerical aspect, 
it would not affect the question posed here, namely: can the LIT be inverted when the response function arises from 
a particle production process? A problem is the slow fall-off of the response function for high W values. From the 
definition of the transform it is clear that the response must behave asymptotically like W 1 ~ x where x > in order 
for the integral to converge. In the case of nuclear photoabsorbtion without pion production the response function 
falls off very rapidly for increasing W resulting in L(gr, gj) also falling rapidly for large gr. Inversion then gives an 
accurate account of 1Z(W) over its entire range by calculating L(gr, gj) for a finite number of gr values. That the 
response function for inclusive photoproduction may have a significantly different asymptotic behavior can be seen 
from the basic Kroll-Rudermann [l3T | cross-section for p(j, n + )n 
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where En are the nucleon energies, E cm is the cm. energy k + E p (k) and a the fine structure constant. This cross 
section approaches a constant for large E cm and the response lZ(E cm ) rises linearly in E cia . In the next section it 
will be seen that the pion production response function for a finite nucleus is tempered at large E cm by structure 
effects. Nevertheless it still rises over a large energy region thereby requiring a different approach for the inversion of 
the transform. 



II. RESULTS 



Details of the model and the parameters used are given in the Appendix. There it is seen that the vector 
I ' l l ) M d ,\is T ii 1 0"/, k p ) ) can also be labeled by the final neutron-neutron spin and that the response function is therefore 
a sum of singlet and triplet contributions. The separable potential used here only has scattering in the spin S = 0, 
isospin T = 1 channel while the S = 1, T = 1 final state consisting of odd partial waves is non- interacting. Fig. 1 
shows the transform Lk p {GR, gj) for the case 07 = k p = 10 MeV. One notes that reaches a maximum at gr ps 200 
MeV and then falls very slowly with increasing gr. Also shown in this figure is the response function TZk (W) for 
kp = 10 MeV. The response function behaves similarly to L as one would expect since L(gr,gi) samples 1Z mainly 
from a region centered at W = gr. We note that lZk p (W) does not rise linearly with energy as in the case of ^(7, Tr + )n 
but falls off slowly with energy after reaching a maximum. Unfortunately, the fall-off of L with gr is so slow and 
covers such a large energy range that inverting it to obtain 1Z(W) over the full range, as was possible in the earlier 
photoabsorption calculations, is not only very difficult but largely physically meaningless because our model is only 
valid for non-relativistic mesons. In fact the present model is only sensible for energies W not exceeding approximately 
40 MeV which corresponds to the maximal pion energy still being non-relativistic. Our aim then is to calculate L 
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only for the segment < or < 40 MeV and then invert it to extract lZk p (W) for W in the same energy range. Since 
TZ can be calculated directly the error can be easily assessed. 

The inversion process we use has already been described in Q but a brief account is as follows. The response 
function TZk p (W) is written as a sum 



N N 



(w) = * w s+2 e- ftw = E n k ( w ) ( 17 ) 



1=1 



and the parameters and $ are determined by fitting 



" r°° (W) 



to L as computed from Eqs. (|12[l - i|13[l of the Appendix. Note that the threshold behavior W s+2 , where S is the 
total spin, is appropriate to the photoproduction process |ll| and differs slightly from the form used in the earlier 
photoabsorption calculations. It turns out however that the quality of fit is nearly independent of which threshold 
behaviour is used. Fig. 2 shows the quality of the fit obtained using N = 10. With the parameters so determined one 
can calculate the response functions and finally the inclusive cross section. As mentioned earlier this model problem 
is simple enough that the response functions are easily calculated in the traditional manner 0, [n| . Fig. 3 shows 
the relative error between the response function as calculated from the LIT technique ( IZat ) and the one calculated 
in the traditional manner. Over most of the energy range between < W < 40 MeV the relative error is in the 1% 
range. However near the threshold, i.e. W < 2 MeV the response function tends to zero for W — ► and limited 
numerical accuracy produces an exaggerated relative error in this region. These effects are too small to be visible in 
a plot of the cross section however. Finally our computed cross section along with data of Booth et al. fijf is shown 
in Fig. 4. The results of using the LIT technique are indistinguishable from the traditional methods in |lOLIll| . 

It is instructive to see the differences in the response functions that occur if one uses transforms calculated in 
different ranges of or. One expects that because of the nature of the Lorentz transform that one should only have 
to calculate the transform up to or = 40 MeV, as was done above, if the response function was only required in 
the energy range < W < 40 MeV. To show this we calculate the transform in four ranges < or < <Ti?, max where 
CTR max = 20,40,60, and 70 MeV. The respective response functions obtained by inverting each of these transforms is 
denoted by lZk p (W|t7fl max ). Fig. 5 shows the error in the erfi max = 20,40,60 MeV cases relative to the cr Rmax = 70 MeV 
case. One notes that the o\R max = 20 MeV case only is accurate for energies up to 20 MeV, that the <JR max = 40 MeV 
case has a less than 1% error at W= 40 MeV, and that the crR max = 60 MeV case has only a very small error up to 60 
MeV. 



III. CONCLUSIONS 



We have shown that the LIT technique can be used to calculate inclusive meson photoproduction cross sections 
in the non-relativistic regime. Rather than fitting the transform over its entire range as was possible in earlier 
photoabsorbtion calculations one fits here only the low energy segment to obtain the low energy response functions. 
The model calculation used here shows that the response functions thus obtained are as accurate as numerical 
techniques will allow. Our next step will be to add a pion-nucleon interaction to the Hamiltonian in order to take 
account of pion scattering effects. The dynamical model of Darwish, Arenhoevel, and Schwamb would provide a 
benchmark calculation against which to further check our method. It would also be of interest to apply these types 
of calculations to higher A nuclei such as 3 He, 3 H, or 4 He. The LIT technique extends readily to these cases as 
well as being able to handle realistic potentials including Coulomb effects. One should expect to be able to study 
meson-rescattering effects with realistic nuclear models in a variety of light nuclei. 
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Appendix 

For the NN interaction we use the separable model of Y. Yamaguchi [12j , which is a pure 5- wave interaction 



V(p,p') = -\)9o(jp)9o(p') j (1 - <?i • <x 2 ) 



9 s(f) 



-Ai5 , i(p)gi(p")^( 3 + (? i -ct 2 ) , 



P 2 + P S 



(1) 
(2) 



Here the labels and 1 refer to the spin-singlet and spin-triplet parts respectively. The following more up to date 
constants for this model have been taken from [T(ij ]: 



a = 0.2316 fin" 1 , ft = 1.129 fm" 1 , ft = 1.392 fni" 1 , 
A = 0.02774 fm- , Ai = ^L^±|il! . 



(3) 
(4) 



These constants fit the deuteron binding energy, the experimental values for the singlet and triplet scattering lengths 
and the singlet effective range. 

Using this separable NN interaction one obtains for the deuteron ground state 



V«ft (« + ft) 3 



1 



(p 2 + a 2 ) (p 2 + /3 2 ) 



(5) 



and a binding energy of 2.224 MeV. 

The Lorentz function ^M d .\{&, k v , p x , p y ) can be decomposed into its spin components as 

4>M d ,\(<r, k p , p x , p y ) = ^2 ^Msi^ fc P> Px, Py) 
S=o,i 



(6) 



where Ms — Md + A with Md and A the deuteron and photon polarization, respectively, and the function 
i ) M s ( a , k pi p x , p y ) is a solution of 



where 



m n 2/z 



^Msi^ k p, Px, Py) - ^S9s{Px) C'M s (p y , A) 

S €{0,1} , 



C M s (Py, A ) = / d P9s{p)^ Ms { (J , K,p,p y ) , 



F%{p) = Mp-^) + {-WMp + &) , 



^M s ,M d ,\ 



rn T[ 



S 1 1 A \\ S \ 
-M s A M d ) 1 1 i 1 



(7) 

(8) 
(9) 
(10) 



We take / 2 /4tt=0.078. 

The constant Cf 4s (p y , A) with respect to p x is 



^M s (Py, A ) — m n ;= 'M s ,M d ,\ J 



j3 9s(p)^l(p) 

2 , 2 

p + 7 — m„ cr 



m„ 7T 2 A s 



ft (ft + V7 2 - m„ ?y 



(11) 



with (m n /fi)/2p 2 = 7 2 (p a ). In the case of 5 1 = 1 this constant vanishes. 
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The solution of the Lorentz-equation therefore is 

7S , , - ^ m « "M s ,M d ,A -^f (Px) + m « 9s{Px) (p y , A) 

^ M >, fc p , Px , ^ = ^ A p 2 +72 _ mwg (12) 

and the transform Lk p ((JR, 07) is 

= £ J] (^Msl^Ms) ■ (13) 
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a R or W [MeV] 

FIG. 1: L kp (<jR,ai) (solid line) and 72. fep (W) (dashed line) shown for 07 = k p = 10 MeV. 
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G R [MeV] 

FIG. 2: Relative error of the fit L§* (ctr, 07) compared to the computed transform Lk p (cm, ai). Here we have used 07 = fc p = 10 
MeV for illustration. 
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FIG. 3: Relative error between the response function TZ^ (W) and 1Zk p (W) computed in the traditional method for ai — fc p 
MeV. 
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FIG. 4: Calculated total cross section for D(7,7r + )nn shown together with the data of fl4|| . The cross sections calculated 
either by the LIT technique or the traditional method [ToL ITU are indistinguishable. 
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